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The properties of two-flavored massless Nambu-Jona-Lasinio model in (l+l)-dimensional R 1 x S 1 
spacetime with compactified space coordinate are investigated in the presence of isospin and quark 
number chemical potentials /i. The consideration is performed in the large N c limit, where N c is 
the number of colored quarks. It is shown that at L = oo (L is the length of the circumference S ) 
the pion condensation (PC) phase with zero quark number density is realized at arbitrary nonzero 
Hi and for rather small values of fj,. However, at arbitrary finite values of L the phase portrait of the 
model contains the PC phase with nonzero quark number density (in the case of periodic boundary 
conditions for quark fields). Hence, finite sizes of the system can serve as a factor promoting the 
appearance of the PC phase in quark matter with nonzero baryon densities. In contrast, the phase 
with chiral symmetry breaking may exist only at rather large values of L. 



I. INTRODUCTION 



Recently, great attention has been paid to the investigation of the QCD phase diagram in terms of baryonic as well 
as isotopic (isospin) chemical potentials. The reason is that dense baryonic matter which can appear in heavy-ion 
collision experiments has an evident isospin asymmetry. Moreover, the dense hadronic/quark matter inside compact 
stars is in general also expected to be isotopically asymmetric. To describe the above mentioned realistic situations, 
i.e. when the baryonic density is comparatively low, usually different nonperturbative methods or effective theories 
such as chiral effective Lagrangians and especially Nambu - Jona-Lasinio (NJL) type models [lj are employed. In this 
way the QCD phase diagram including chira l sy mmetry restoration @-@|, color superconductivity (CSC) and 
charged pion condensation (PC) phenomena [10l4l3| were investigated under heavy- ion experimental and/or compact 
star conditions, i.e. in the presence of such external conditions as temperature, chemical potentials and possible 
external (chromo)magnetic fields (see the above references). 

Moreover, it was also realized in the framework of NJL-type models that different finite size factors, i.e. curvature 
or nontrivial spacetime topology as well as a finite spatial volume of a system, can significantly change the properties 
of both the chiral- and color superconductivity phase transitions. This conclusion is the result of chiral symmetry 
breaking investigations in weakly curved spaces [3 UH and in spaces with nontrivial topology, where one or more 
space coordinates are compactified [n| . In addition, the properties of finite size normal quark matter droplets in the 
language of the MIT-bag model were considered, e.g., in the review [l7| ■ It was also noted that the position of the chiral 
critical end point of the QCD phase diagram, which could be investigated in heavy ion collision experiments, depends 
essentially on the size of a finite system [18j . Next, the effect of spaces with constant curvature or nontrivial topology 
on CSC was studied in [l9l |2fij |. In particular, there it was shown that in a restricted region the CSC phenomenon 
might occur much easier than in an infinite one (20| . Moreover, the stability of finite size quark matter droplets in 
the color-flavor locked phase was considered in the framework of a bag model using the so-called multiple expansion 
method [21| . However, up to now we have only scarce information about the formation of the pion condensation 
phase in dense baryonic matter under the influence of finite-size effects (see below) . This fact is partially explained 
by the complexity of the problem arising from the necessity to introduce an additional isotopic chemical potential for 
the adequate description of quark matter with isospin asymmetry and, in particular, of the PC phenomenon. 

Obviously, the (3+l)-dimensional NJL models depend on the cutoff parameter which is typically chosen of the order 
of one GEV, so that the results of their usage are valid only at comparatively low energies, temperatures and densities 
(chemical potentials) . Besides, there exists also a class of renormalizable theories, the (I+l)-dimensional chiral Gross- 
Neveu (GN) type models H3,E1 that can be used as a laboratory for the qualitative simulation of specific properties 
of QCD at arbitrary energies. Renormalizability, asymptotic freedom, as well as the spontaneous breaking of chiral 
symmetry (in vacuum) are the most fundamental inherent features both for QCD and all GN type models. In addition, 
the fi — T phase diagram is qualitatively the same for the QCD and GN model [24h"27| (here \i is the quark number 
chemical potential and T is the temperature). Note also that the GN type models are suitable for the description of 
physics in quasi one-dimensional condensed matter systems like polyacetylene [281 ] . It is currently well understood (see, 
e.g., the discussion in [25U27I ] ) that the usual no-go theorem (29| . which generally forbids the spontaneous breaking 
of any continuous symmetries in two-dimensional spacetime does not work in the limit N c — > 00, where N c is the 
number of colored quarks. This follows from the fact that in the limit of large N c the quantum fluctuations, which 



1 Below wc shall use the notation "NJL2 model" instead of "chiral GN model" for (l+l)-dimensional models with a continuous chiral 
symmetry, since the chiral structure of the Lagrangian is the same as that of the (3+l)-dimensional NJL model. 
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would otherwise destroy a long-range order corresponding to a spontaneous symmetry breaking, are suppressed by 
1/N C factors. Thus, the effects inherent for real dense quark matter, such as color superconductivity (spontaneous 
breaking of the continuous color symmetry) or charged pion condensation (spontaneous breaking of the continuous 
isospin symmetry) might be simulated in terms of a simpler (l+l)-dimensional GN-type model, though only in the 
leading order of the large N c approximation (see, e.g., [26l |3fJ and I3U 1331 . respectively). Finally, one should recall 
that both the chiral phase transition [33[ and color superconductivity |30l| were investigated in the framework of GN 
models with account of the nontrivial R 1 x S 1 spacetime topology. 

In general, this paper is devoted to the consideration of the charged pion condensation phenomenon under the 
influence of finite-size effects. The problem was partially solved in J34J, where PC was studied in the framework of 
the (3+l)-dimensional NJL model in the Einstein universe with a constant curvature and finite spatial volume, and 
in [3l|, where it was considered in the SUl(2) x SUr(2) symmetric GN model in spacetime with nontrivial topology, 
i.e. on the R 1 x S 1 manifold with compactified space coordinate, and in the assumption that N c — > oo. For technical 
simplifications, in those papers the quark chemical potential \i was assumed to be zero, i.e. the considerations were 
performed at nonzero isospin /i/ chemical potential only (this situation corresponds to quark matter with zero baryon 
density). Obviously, the considered problem deserves further, more detailed investigations, this time when both 
chemical potentials are nonvanishing, /J / 0,/jj ^ 0. This is physically motivated by the fact that quark matter, 
which might be created in heavy-ion collisions, has in general nonzero baryon- and isospin densities and must be 
investigated in the framework of a theory with nonzero fj, and Moreover, quark- matter lumps occupy a finite 
volume. Hence, in this paper we study, for illustration, the PC phenomenon in a SUl{2) x SUr(2) symmetric two- 
dimensional NJL2 model with [i ^ 0, fii ^ 0, when the spatial coordinate is compactified and N c — > oo. In particular, 
we shall demonstrate that the finite size of the system promotes the appearance of the PC phenomenon in dense 
baryonic matter. We hope that the results of such a study of a renormalizable two-dimensional model may provide 
an additional stimulus for further investigations of the charged PC phenomenon in more realistic models. 

The paper is organized as follows. In Section II we derive, in the leading order of the large iV c -expansion, the 
expression for the thermodynamic potential of the two-flavored massless NJL2 model with quark number chemical 
potential \i and isospin chemical potential for zero temperature in R 1 x R 1 - and in R 1 x S 1 spacetimes.Then, in 
Section III the phase structure of the model is investigated both at a finite value of L and at L —> 00, where L is the 
radius of the S 1 circumference. Finally, Section IV presents a summary and some concluding remarks. 



II. DESCRIPTION OF THE MODEL AND ITS THERMODYNAMIC POTENTIAL 



A. Effective action 



We consider a (l+l)-dimcnsional NJL 2 model to mimic the phase structure of real dense quark matter with two 
quark flavors (u and d quarks). Its Lagrangian has the form: 



£q.7j — q 



G_ 



[{qqf + (qh 5 rq) 



(1) 



where each quark field q(x) = qi a (x) is a flavor doublet (i — 1,2 or i = u,d) and color iV c -plet (a — 1,...,N C ). 
Moreover, it is a two-component Dirac spinor (in ([IJ the summation over flavor, color, and spinor indices is implied); 
Tk(k = 1,2,3) are Pauli matrices. The quark number chemical potential \i in ([T]) is responsible for the nonzero 
baryonic density of quark matter, whereas the isospin chemical potential /i/ is taken into account in order to study 
properties of quark matter at nonzero isospin densities (in this case the densities of u and d quarks are different). The 
Dirac gamma matrices in two-dimensional spacetime have the following form: 
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(2) 



Let us consider the symmetries of the Lagrangian. If jij = 0, the Lagrangian ([T]) is not only SU(N C ) symmetric, but 
also invariant under transformations of the chiral SUl(2) x SUr(2) group. However, if fii ^ 0, the latter symmetry is 
reduced to Ui 3 l{1) x Ui 3 r(1), where ^3 = T3/2 is the third component of the isospin operator (here and above the sub- 
scripts L, R mean that the corresponding groups only act on left- and right-handed spinors, respectively). Obviously, 
this symmetry can also be presented as Ui a (l) x [7^/ 3 (l), where Z7/ 3 (l) is a vector isospin subgroup, and Uai 3 0-) is 
an axial isospin subgroup. Quarks are transformed under these subgroups as q — > exp(zar3)g and q exp(iaj 5 T3)q, 
respectively. 

The linearized version of Lagrangian ([TJ, which contains composite bosonic fields o{x) and 7r a (a;)(a = 1, 2, 3), has 
the following form: 



£<7,7r — q 



7 ly i<9„ + ^7° + tt t 37° - c - i7 5 7r Q T a 



Nc r m 1 
q ~ 4G ^ 7r ° 7ra ^ 



(3) 
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^From Lagrangian Q one obtains the following constraint equations for the bosonic fields: 

G G 
<*{x) = - 2 ]v^(^); Ta(ar) = -2^-(gi 7 5 r Q g). (4) 

Obviously, Lagrangian ([3]) is equivalent to Lagrangian ([1]), when using the constraint equations ((!]). Furthermore, it 
is clear from ((U) that the bosonic fields are transformed under the isospin Ui 3 (2) and axial isospin Uai 3 (2) subgroups 
in the following manner: 

Uj s (l) : o — y o\ 7T3 — > 7T3; 7ri — ► cos(2a)7ri + sin(2a)7T2; 7r2 — > cos(2a)7T2 — sin(2a)7ri, 
tr,4j 3 (l) : 7i"i — > 7Ti; 7T2 — > 7T2; cr — > cos(2a)cr + sin(2a)7T3; 7r3 — > cos(2a)7T3 — sin(2a)<7. 

To avoid the no-go theorem, which forbids the spontaneous breaking of continuous symmetries in the considered 
case of one space dimension, we restrict the discussion only to the leading order of the large N c expansion (i.e. to the 
case N c — > oo), where this theorem is not valid (25l - [27| . In particular, the effective action S e g[a, n a ] can be found in 
this approximation through the relation: 



= f[dq][dqV Id2xL ^. (6) 

Carrying out the Gaussian-type path integration in (|6]) over fermion fields, we obtain 

e iS«tt[<T,,r a ] = e -i£f /rf 2 ^ 2 +O detD . (7) 
In ([7]) we have used the notation D = D x 7 C , where I c is the unit operator in the iV c -dimensional color space and 

D = vfdu + /i 7 ° + ^-T 37 - a - i 7 5 ^aT a . (8) 

Then using the general formula detD = exp Tr cs /a;lnD and taking the logarithm of ([7]), one obtains the following 
expression for the effective action in the leading order of the 1 /./Vc-expansion: 

(•2,2 

S e s{<7, 7T a ) = —N c / 4G a d 2 x - iN c Tr sfx ]nD, (9) 

where the Tr-operator stands for the trace in spinor (s), flavor (f), and two-dimensional coordinate (x) spaces, 
respectively (the trace in color (c) space, obviously, equals N c ). 



B. Thermodynamic potential 

Starting from ((9|), one can define the thermodynamic potential (TDP) of the model at zero temperature T = in 
the mean- field approximation, i.e. in the leading order of the l/iV c -expansion |3lj : 



O (rr «r \ ~ ^(o", 7T a ) 



N c Jd 2 



x 



4G f d 2 
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ft 



Tr s/x lnD 
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where the mean fields a and 7r a are x— independent quantities, and in the round brackets of (|I0[) just the momentum 
space representation D = jp + /i 7 ° + ^t^j — a — i 7 5 7r a T Q of the Dirac operator D appears. □ Obviously, Tr s ^lnD = 
Y) . lne^ , where summation over all four eigenvalues e$ of the 4x4 matrix D is implied and 



ei.2,3,4 = -a± \l(p + fJ,) 2 -p 2 -n 2 +[^Y± ^ (p + fi) 2 - tt? - Trf. ( ' ' I 

Hence, 



d 2 p 

AG 'J (2tt) 2 



S2 MMJ (cr,7r a ) = — — hi / 7— 2ln(eie 2 e3e4j 



(12) 



2 There exists also the interesting case of x-dependent mean fields describing chiral density waves [35ll , For the following investigation of 
finite-size effects we restrict us, however, for simplicity, to the case of constant fields. 
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where 




(13) 



It is clear that the TDP J7 MMJ (tr, 7r a ) is symmetric under the transformations [i — > — /i and/or [ii — > —fii. So it is 
sufficient to consider only the region fi > 0, fj,i > 0. Taking into account this constraint and integrating in (|12p . one 
obtains the following expression for the TDP of the system: 

2 2 °° 

n M wfa*-a) = a ~ J ^{e+ + £- + (/i -e+)0(ji-e+) + (ji - £_)0(/x - £_)}, (14) 

— oo 

where 0(x) is the Heaviside step function. To simplify the task, let us note that both the quasiparticle energies (fTTj) 
and hence the TDP (fT4|) depend effectively only on two combinations (w 2 + 7r|) and (7r| + a 2 ) of the bosonic fields, 
which are invariants with respect to the Ui 3 (l) x {7at 3 (1) group, as is easily seen from ([5]). In this case, without loss 
of generality, one can put 7T2 = = 0, and study the thermodynamic potential (| 14[) TDP as a function of only two 
variables, M = a and A = tti. Then the TDP has the following form: 

oo 

M 2 + A 2 f dm 

n m (M,A) = ^^ J -^{E+ + E A + (f,-E+)9^-EX) + ^-E A )9^-E A )}, (15) 

— OO 

where 



E£ = ^[\/Pl + M 2 ±f) + A^^{E±v) 2 + A 2 - \E = ^p\ + M 2 , „ = f J . (16) 
Since we are going to study the phase diagram of the initial NJL2 model, the system of gap equations is needed: 

<9M ' 9A V ; 

The coordinates M and A of the global minimum point (GMP) of the TDP (fl"5|) provide two order parameters (gaps), 
which are proportional to the ground state expectation values < qq > and < qi^T\q >, respectively. In this case the 
gap M is just a dynamical quark mass. Obviously, the pair of gaps M and A is a solution of the system (|17l) . So, the 
GMP of the form (M = 0, A = 0) corresponds to the phase, in which the initial J7j 3 (l) X Uaj 3 (1) symmetry remains 
intact. If M ^ 0, A = in the GMP, then the Z7V 3 ( 1 ) symmetric phase is realized in the model. Finally, the GMP 
of the form (M = 0,A ^ 0) corresponds to the Uai 3 (1) symmetric phase, where charged pions are condensed and 
isospin symmetry, C7/ 3 (l), is broken spontaneously. (Note that due to the zero value of the bare quark mass in the 
Lagrangian (1), the GMP of the form (M ^0,A/0) does not appear for the TDP (fTS")) .) On the basis of the gaps 
M and A the following quantities can be introduced, 

Pou = E A -fi, p od = E^-iJ,, pm = -(E A ~+n) p ol = -(E A + fx), (18) 

which are the energies of u-, d-, u-, d-quasiparticles (quarks) or one-fermion excitations of the corresponding phase. 
The quantities po u ,Pod from (|18|) are the energies necessary for creation of u and d quarks with momentum pi, whereas 
Pott,Pod are t ne en ergies necessary for annihilation of u and d antiquarks. 

It is clear that the TDP (|15l) is an ultraviolet divergent quantity, so one should renormalize it, using a special 
dependence of the bare quantities, such as the bare coupling constant G = G(A), on the cutoff parameter A (A 
restricts the integration region in the divergent integrals, |pi| < A). The detailed discussion of this procedure was 
performed, for example, in papers [3lT |32|] . The main stages of it are the following. First, we transform the expression 
for the TDP OB, 



JWM, A) = V (M, A) - J ^{Ei +E A - 2^J p\ + M 2 + A 2 + (/i - E£)8(^ - E+) + ( M - E A )9( t , - E A )}. 

— OO 

(19) 

The integral in the above expression is convergent and Vq(M, A) looks like 



%(M,A) = fi MMJ (M,A) = - + A -2 / p-Jtf + AP + A*. (20) 
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Second, it is evident that the effective potential ([20]) is a divergent quantity. So, to renormalize it we cut the 
integration region in (f2"U|) . \px\ < A, and require the bare coupling constant G = G(A) to have the following form 

xm = *h dpi 7Wrwr« { — ^ — )■ (21) 

Then in the limit A — > oo one can obtain a finite, i.e. renormalized, expression for Vq(M, A): 



, . M 2 + A 2 r 
Vb(M,A = — 5- 

Z7T 



. M' , _ . 



(22) 

J o / 

Since Mo might be considered as a free model parameter, it follows from (|2Tj) and (|22[) that the renormalization 
procedure of the NJL2 model is accompanied by the dimensional transmutation phenomenon. Indeed, in the initial 
unrenormalized expression (j2"0")l for Vo(M, A) the dimensionless coupling constant G is present, whereas after renor- 
malization the effective potential (|2"2"j) is characterized by a dimensional free model parameter Mq. Moreover, as it 
is clear from (|22|) . the global minimum point of the affective potential Vo(M, A) lies just at the point M = Mq. So 
in vacuum, the chiral SUl(2) X SUr(2) symmetry of the NJL2 model (1) is always spontaneously broken and the 
quantity Mq might be treated as a dynamical quark mass (in vacuum). 

Now, taking into account the expression (|2"2")l for Vq(M, A), we see that the TDP (Tl"9")) is a finite renormalization 
invariant quantity which describes the properties of dense and isotopically asymmetric quark matter in an infinite 
volume. 

For the further analysis we need also the expressions for the quark number density n q and isospin density ni in the 
phase with gaps M and A which follow directly from (|19p . 

00 



(23) 



— OO 

OO 



C. Compact ificat ion of the spatial coordinate 



Since physical effects generally take place in a restricted space region, in the present paper we are going to study 
the influence of a finite size of the system (at zero temperature) on the thermodynamical properties of dense and 
isotopically asymmetric quark matter. To simulate a real situation, we put our (l+l)-dimensional system with 
Lagrangian ^ into a restricted space region of the form < x < L (here x is the space coordinate). Then a 
corresponding modification of the TDP ([T5]) is needed. 

It is well-known that such a constraint on any physical system is equivalent to its investigation in a spacetime with 
nontrivial topology in which space coordinates are compactified. In our case it means that we can consider the model 
(1) in spacetime with the topology R 1 x S 1 and with quantum fields satisfying some boundary conditions of the form 

q(t,x + L) = e m q(t,x), (24) 

where < a < 2, L is the length of the circumference S , and the variable x means the path along it. Below, we shall 
use only two values of the parameter a : a — for a periodic boundary condition and a = 1 for the antiperiodic one. 

As a consequence, to obtain the TDP Flium (M, A) of the initial system placed in the restricted domain < x < L 
and at nonzero chemical potentials /i 7^ 0, /i/ 7^ 0, one must simply replace the integration in (|19[) and (1201) by an 
infinite series, according to the rule: 



00 

/d 1 
^/(Pi)->£ E /(Pm)i Pm = £( 2 " + «). n = 0,±l,±2,. 



(25) 



As a result, we obtain for the corresponding TDP ^^(Af, A) in the spacetime of the form R 1 x S 1 the following 
expression 



1 — / 2 

n=-oo (26) 

+ (P - Ei An )6((i - E+ An ) + (p- E- LAn )6{p - EZ An )}, 



3 Recall, the thermodynamic potential in vacuum, i.e. at /Lt = fXi = 0, is usually called effective potential. 
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where p = \/M 2 + A 2 



and 
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LAn 



\ 



' M 2 + j^(2n + a) 2 ± 1/ 



A 2 



VL(p)-Vl(0) = 
2 °° 



V 47T 



7T 



P 2 7 



2 

- — jVp 2 L 2 + n 2 a 2 



2na 



y / Ti 2 (2n + a) 2 + L 2 p 2 + v / 7r 2 (2n - a) 2 + L 2 p 2 - Ann 



p 2 L 2 
2mr 



(27) 



(28) 



Here 7 = 0.577... is the Euler constant. The expression (|28fl is the generalization of the vacuum effective potential 
([2"12]) in the spacetime with trivial R 1 x i? 1 -topology to the case of s pac etime with nontrivial topology of the form 
R 1 x S 1 (the detailed derivation of (|28|) is presented in Appendix of [31|). For further discussion, we need also the 
expression for the quark number density n q L in the R 1 x S 1 spacetime, which can be easily obtained from the TDP 



(29) 



Moreover, in what follows it will be convenient to use the dimensionless quantities: 



A = 



LMo 



JL 

Mo' 



v _ Hi 
M ~ 2Af ' 



M 
Mo' 



S = 



M ' 



(30) 



III. PHASE STRUCTURE OF THE MODEL 



In the present section we shall study the phase structure of the NJL2 model (1) depending on three external 
parameters, p, v and A (|30l) in the case T = 0. In general, this first means that for each fixed values of ft, v and A we 
should study the global minimum point of the TDP Qx^j (M, A) as a function of M, A (or dimensionless variables 
to, S (|30[) ) and then indicate the symmetry group of this point. Secondly, we must divide the three-dimensional space 
(m > 0, v > 0, A > 0) into regions (phases), where each of them is composed by the points corresponding to the 
same symmetry group of the TDP global minimum point. Since this is a rather difficult problem, one can use an 
equivalent way. Instead of studying the above mentioned three-dimensional phase diagram, we will slice it into the 
planes A = const and then perform the above procedure only in the planes labeled by some values of A. (Of course, 
the more values of A are taken into consideration, the more exact is the obtained picture about the phase structure 
of the model.) 



A. Particular case: R 1 x R 1 spacetime (A = 0) 

Let us begin our study of the (z>, fl, A)-phase diagram of the model (1) starting from the plane A = (L = 00) 
that corresponds to the spacetime of trivial topology R 1 x R 1 . In this case the numerical investigations of the global 
minimum point (GMP) properties of the TDP (|T9")) bring us to the following phase structure of the model shown in 
Fig. 1. There, in the corresponding (D, /i)-plane, one can see four phase regions denoted by PC as well as by the 
roman figures I, II and III, respectively. 

In the region denoted by PC, i.e. in the charged pion condensation phase, each point (p, p) corresponds to the GMP 
of the form (to = 0, i5 = 1). Hence, everywhere in this phase the relation jl < fi c < 5 is valid, since p c = l/y/2. As a 
result, we see from (| L8[) that all quasiparticles are gapped in the PC phase, i.e. a finite amount of energy is needed 
to create both u- and d-quarks. Moreover, it is easily seen from (|23p that the quark number density n q is identically 
zero in this phase. In contrast, the isospin density ni is nonzero and equal to v /it at each of the PC phase points. 

Then, let us observe the properties of phases denoted by I, II, III. At each point of these phases the gap 5 — 0, but 
the gap to = appears only in the phases I and III. As a result, in regions I and III the chirally J7/ 3 (l) X Uai 3 0-)- 
symmetric phase with massless quarks is arranged. Contrary, in the region II, where the order parameter to is nonzero, 
this symmetry is spontaneously broken down to the isospin Ui 3 subgroup. In this region the order parameter to is a 
smooth function vs p, and v and tends to zero when p and v tend to infinity (to — > when p — ¥ 00 and v — > 00). We 
call this phase the normal quark matter phase, since here quarks dynamically acquire a mass which is equal to the 
order parameter to (multiplied by Mo), and space parity is not broken. Note also that the upper and lower boundaries 
of the region II tend asymptotically to the line p = v from the upper- and lower sides, respectively. In the regions I 
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I 










III 










PC 

1 1 




Symmetric phase 



PC 



FIG. 1. The phase portrait of the model at 

A = and v > 0. In the phases I and III the initial 
Ui 3 (l) x Uai 3 {1) symmetry is not broken. In the phase 
II it is broken down up to Ui 3 (l), and in the pion con- 
densation (PC) phase with zero quark number density it 
is spontaneously broken down to Uai 3 (V). 



FIG. 2. The periodic case: The (y, \x) phase portrait at 
A = 0.1. PCd and PC mean the pion condensation phases 
with nonzero and zero quark number densities, respec- 
tively. II is the chirally non-symmetrical phase. 
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Symmetric phase 



PCd L 




Symmetric phase 



PCd 



PC 
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FIG. 3. The periodic case: The (v, fi) phase portrait at 
A = 1. PCd and PC mean the pion condensation phases 
with nonzero and zero quark number densities, respec- 
tively. 



FIG. 4. The periodic case: The {v, p) phase portrait at 
A = 2. PCd and PC mean the pion condensation phases 
with nonzero and zero quark number densities, respec- 
tively. 



and III, as it follows from (|23[) . the quark number density is n q = 2/z/7r and the isospin density is nj = is/ir. Finally, 
we would like to discuss the difference between I and III phases. It is evident that in phases I and III the dispersion 
relations (|18[) for quasiparticles have the following simple form: 



Pou = \pi -v\- n; pa cl = Pi + v - jj,. 



(31) 



Therefore, it is easy to check that in phase I both u- and d- quasiparticles are gapless. This means that to create these 
quarks costs no energy, i.e. at each fixed values of /i, v there exist space momenta p\ and p\* such that Pou(p*) = 
and podip**) = 0. In contrast, in phase III only it-quasiparticles are gapless, but d-quarks are gapped. 
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B. General case: R 1 x S 1 spacetime (periodic boundary conditions) 

In the present section we consider both the phase structure and relevant properties of the model (1) at some 
particular nonzero values of A in the case of periodic boundary conditions, i.e. at a = in (1241) . In terms of the 
dimensionless quantities (|30p the TDP (j2"o| has in the periodic case the following form: 

4o iW! (M, A) = (m 2 + <5 2 ) [ln(4A) - 7 ] - A^ (m + £) 2 + <5 2 - Ay 7 (to — i>) 2 + <5 2 

" 2A E { £ n«=o + ^wo ~ 4nA - I?L±- + (m- - S+ a=0 ) + (£- En*=o)0fr - E~ a= o)}, (32) 

n=l 

where 7 = 0.577... is again the Euler constant and 

£4=0 = \/{y/m 2 + (2nXf±P) 2 +S 2 . (33) 

We also need an expression for the quark number density n q L in the periodic case which follows directly from (|29[) at 
a = 0, 

h qL = -^-;-n qL = \9(fi - y/(m + i>) 2 + S 2 ) + \6(fi - y/(m - i>) 2 + 5 2 ) 
Mq 

00 

+ 2XJ2{0^ ~ E+ a=0 ) + 6{fi - E- a=Q )}, (34) 

71=1 

where the zero mode terms are selected in an explicit form. The results of numerical investigations of the global 
minimum point (GMP) of the TDP (f3"2"j) are presented in Figs. 2-4, where the phase portraits of the model are 
depicted at A = 0.1, 1, 2, respectively. There, in the symmetric phase the global minimum of the TDP lies at the 
Uj 3 (l) x Uaj 3 (1) symmetric point (m = 0,(5 = 0). The phase II corresponds to the Uj 3 (l) symmetric GMP of the 
form (to ^ 0, 5 = 0). In this phase quarks are massive, and the quark number density n q L is nonzero. In contrast 
to the case A = (L = oo), we see that at A ^ the phase II occupies a compact region in the phase diagram and 
completely vanishes at A > A p ~ 0.16. 

Moreover, one can see in Figs 2-4 two pion condensation phases, PC and PCd. They correpond to the GMP of the 
TDP (|3"2"j) with (m = 0,(5 7^ 0) and the Uai 3 (X) symmetry group. Clearly, the gap 5 of the pion condensation phases, 
as well as the gap to in the above mentioned phase II, depend on the external parameters, i.e. 5 = 5(p,, v 1 A) (see 
Figs 5, 6). The main difference between these pion condensation phases is the following. In the PC phase the quark 
number density n q ^ is equal to zero, whereas in the PCd phase this quantity is nonzero. The boundary between PC 
and PCd phases is defined by the zero mode term in (|34[) at m = and 5 = S(fb, v, A), i.e. it is the plot of the function 
jl = fi{p, A) which is defined implicitly by the equation 

A = y/i> 2 + 6 2 (fi,D,X). (35) 

Note, we have found the PCd phase at each arbitrary small A > 0, whereas at A = 0, i.e. in the ordinary (1+1)- 
dimensional Minkowsky spacetime, it is absent. This means that a nontrivial spacetime topology promotes the creation 
of the charged pion condensed phase with a nonzero quark number density. The behavior of both the gap 8{fi, A) 
and quark number density n q L (|34[) vs jl in the case of a periodic boundary condition is illustrated in Fig. 5 at A = 2 
and v = 0.1. 



C. General case: R 1 x S 1 spacetime (antiperiodic boundary conditions) 

In the antiperiodic case, i.e. at a = 1, the TDP (|26|) has the following form in terms of variables (|30|) : 

00 2 + (5 2 

j£jn iwiI (M, A) = (m 2 + <5 2 )[ln(A) - 7 ] - 2A £ {e+ =1 + E~ =1 - 2(2n + 1)A - — 

n=0 ' (36) 

+ (A - K=xW - K=x) + (A - E- =1 )6{p, - E- =1 )}, 

where 

B± =1 = yj (v/m 2 + (2n + 1) 2 A 2 ± v) 2 ~+7 2 . (37) 

The investigation of the TDP (|36l) leads to the following conclusions about the phase structure of the model (1) 
(see Figs 7-9): First, as in the case with periodic boundary conditions for quark fields, at A ^ oo the chirally non- 
symmetrical phase II in the antiperiodic case occupies a compact region in the (v, /i)~plane. In addition, this phase 
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FIG. 5. The periodic case: The behavior of the gap 8 
(curve 1) and quark number density h q L (|34|l (curve 2) 
vs ft at A = 2, v = 0.1. 



1.4 



0.6 





FIG. 6. The periodic case: The behavior of the gap 5 vs 
v at A = 2 and arbitrary /} £ (0, 0.7). 




FIG. 7. TVie antiperiodic case: The phase portrait at 
A = 0.1. PC again means the pion condensation phase 
with zero quark number densities. II is the chirally non- 
symmetric phase. 



FIG. 8. The antiperiodic case: The phase portrait at 
A = 1. PC means the pion condensation phase with zero 
quark number densities. 



is realized only at A < A a m 0.66. Secondly, in the antiperiodic case the phase portrait of the model contains only 
the pion condensed phase with zero quark number density at arbitrary A-values, i.e. a PC phase with nonzero quark 
number densities is absent. The behavior of the gap 6 vs v at A = 2 and arbitrary fixed value of ft £ (0, 1) is depicted 
in Fig. 10. 

Finally, note that the above results, obtained both in the periodic and antiperiodic cases, refer strictly to the case 
v > 0. A detailed study of GN-type models in spacetime with nontrivial topology in the case of [i =t 0, v — was 
made in (33j. 



IV. SUMMARY AND DISCUSSION 



In this paper we have studied the so-called charged pion condensation (PC) phenomenon which might be observed 
in dense baryonic matter. Since we were going to pay special attention to the role of finite volume and to concentrate 
our attention on the principle properties of this effect, we have restricted ourselves to the (l+l)-dimensional NJL2 
model (1) with finite values of quark number /1 and isospin /j,t chemical potentials at zero temperature. Moreover, 
our consideration was performed in the leading order of the large iV c -expansion. 
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1.4 



Symmetric phase 




FIG. 9. The antiperiodic case: The phase portrait at 
A = 2. PC means the pion condensation phase with zero 
quark number densities. 




FIG. 10. The antiperiodic case: The behavior of the gap 
8 vs v at A = 2 and arbitrary jl G (0, 0.7). 



Recall that the charged PC phenomenon was studied recently in the framework of some QCD-like effective theories 
such as NJL-type models or chiral effective theories in usual (3+l)-dimensional Minkowski spacetime [pol - tl3| . However, 
the existence of the PC phase with nonzero baryon (quark number) density, which we have called the PCd phase, was 
predicted there without any confidence. Indeed, for some values of model parameters (the coupling constant G and 
cutoff parameter A, etc.) the PCd phase is allowed by the NJL-type models. However, for other physically interesting 
values of G and A the PCd phase is forbidden in the framework of these models [HI]. Moreover, if the electric charge 
neutrality constraint is imposed, the pion condensation phenomenon depends strongly on the bare (current) quark 
mass values. In particular, it turns out that the PCd phase is forbidden in the framework of the NJL-type models if 
bare quark masses reach the physically acceptable values of 5 -j- 10 MeV (see [Hj]). In addition, recent investigations 
of the charged pion condensation phenomenon in terms of the (l+l)-dimensional massive/massless NJL 2 model in 
R 1 x R 1 spacetime shows that the PCd phase is also absent there [32| (see also Sec. IIII Al of the present paper). 

In this paper we have placed the NJL 2 system (1) into a restricted spatial region or, equivalently, in the spacetime 
with nontrivial topology of the form R 1 x S 1 , where the space coordinate is compactified into a circumference of a finite 
length L. It turns out that at A = (A ~ 1/L) the (jj,t, /i)-phase diagram of the model contains the zero quark number 
density PC phase as well as two symmetric phases, where the initial symmetry of the model remains intact, and the 
phase II with broken chiral symmetry (see Fig. 1). Then, if A > 0, we have observed in the (/xj, /i)-phase diagram a 
more physically interesting situation, since at periodic boundary conditions for quark fields and for rather small values 
of [ii the PCd phase is realized there (see Figs 2-4). So we may conclude that finite size promotes the appearence of 
the PC phase with nonzero baryonic density at least in the framework of the NJL 2 model (1). This is one of the main 
results of our paper and we hope that such an effect also takes place in the more realistic case of (3+l)-dimensional 
spacetime. Taking into account the fact that Cooper pairing of quarks is also significantly facilitated at L ^ oo [2Cj| . 
it is reasonable to expect that finite size can affect key properties of any system in comparison with those obtained 
from the L — > oo limit (see, e.g., the papers [13, EH, E( , where this point of view is also supported). 

Among other results of our paper, it is interesting to note that the chirally broken phase II exists on the (fij,fi)- 
phase diagrams only at A < X p ss 0.16 (A < A a « 0.66) in the case of periodic (antiperiodic) boundary conditions. In 
contrast, the PC phase (PCd phase in the case of periodic boundary conditions) is realized at an arbitrary value of 
A > 0. 

Obviously, it would be also interesting to study the influence of finite-size effects on the existence of chiral density 
waves and on the PC phenomenon in some QCD-like models in (3-(-l)-dimensional spacetime. 
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